Low temperature shape relaxation of 2-d islands by edge diffusion 
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We present a precise microscopic description of the limiting step for low temperature shape relaxation 
of two dimensional islands in which activated diffusion of particles along the boundary is the only 
mechanism of transport allowed. In particular, we are able to explain why the system is driven 
irreversibly towards equilibrium. Based on this description, we present a scheme for calculating 
the duration of the limiting step at each stage of the relaxation process. Finally, we calculate 
numerically the total relaxation time as predicted by our results and compare it with simulations of 
the relaxation process. 
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I. INTRODUCTION 

The understanding, description and control of struc- 
tures at the nanometer scales is a subject ofrinterest from 
the fundamental and applied points of viewuH. From 
fundamental point of view, there is a large literature 
concerning the growth of crystals and their shape. Yet, 
while the description of the equilibrium shape is rather 
clear, the dynamic description of crystal growth is still 
not well understood. In particular, we lack a complete 
understanding of the time scales involved in the relax- 
ation process, and the mechanisms which irreversibly 
conduce the island to its equilibrium shape. 

In this work, we study the shape relaxation of two di- 
mensional islands by boundary diffusion at low tempera- 
tures. The typical size of the islands we will be concerned 
with consists of a few thousand atoms or molecules, cor- 
responding to islands of a few nanometers. The model 
we consider is the same as the one studied irjj, where 
two mechanisms of relaxation, depending on tempera- 
ture, were pointed out: At high temperatures, the clasa 
sical theory developed by Herring, MuUins and NicholsS 
appears to describe adequately the relaxation process. 
In particular, it predicts that the relaxation time scales 
as the number of atoms to the power 2. However, at 
low temperatures, the islands spend long times in fully 
faceted configurations, suggesting that the limiting step 
of the relaxation in this situation is the nucleation of a 
new row on a facet. This assumption leads to the cor- 
rect scaling behavior of the relaxation time on the size 
of the island, as well as the correct temperature depen- 
dence. Yet, it is unclear what drives the island towards 
equilibrium in this scenario. 

In this paper we propose a detailed description of 
this low temperature relaxation mechanism, and identify 
the event that drives the island towards its equilibrium 
shape. Based on our description, we construct a Markov 
process from which we can estimate the duration of each 
stage of the relaxation process. Finally, we use our re- 



sult to determine the relaxation time of the islands and 
compare with simulation results. 

The specific model under consideration consists of 2D 
islands having a perfect triangular crystalline structure. 
A very simple energy landscape for activated atomic mo- 
tion was chosen: the aim being to point out the basic of 
mechanisms of relaxation, and not to fit the specific be- 
havior of a particular material. The potential energy Ep 
of an atom is assumed to be proportional to the number 
i of neighbors, and the kinetic barrier Eact for diffusion 
is also proportional to the number of initial neighbors 
before the jump, regardless of the final number of neigh- 
bors: Eact — —Ep = i * E where E sets the energy scale 
(_E = 0.1 cV throughout the paper). Therefore, the prob- 
ability Pi per unit time that an atom with i neighbors 
moves is pi = z/oexp[— i * E/ki,T], where vq — lO^^s^^ is 
the Debye frequency, ks is the Boltzmann constant and 
T the absolute temperature. Hence, the average time in 
which a particle with i neighbors would move is given by 



- exp[i * E/hT] 



(1) 



The complete description of the model and of the simu- 
lation algorithm can be found iro, where it was studied 
using Standard Kinetic Monte Carlo simulations. This 
simple kinetic model has only one parameter, the ra- 
tio E/ksT. The temperature was varied from SS/iT to 
hQQK , and the number of atoms in the islands from 90 
up to 20000. The initial configurations of the islands 
were elongated (same initial aspect ratio of about 10), 
and the simulations were stopped when the islands were 
close to equilibrium, with an aspect ratio of 1.2. The 
time required for this to happen was defined as the re- 
laxation time corresponding to that island size and tem- 
perature. Concerning the dependence of the relaxation 
time on the size of the island, two different—behaviors 
depending on temperature were distinguishecO. At high 
temperature, the relaxation time scaled as the number of 
atoms to the power 2, but this exponent decreased when 



temperature was decreased. A careful analysis showed 
that the exponent tends towards 1 at low temperature. 
The dependence of the relaxation time as a function of 
temperature also changes, the activation energy was 0.3 
eV at high temperature and 0.4 eV at low temperature. 
In this context, it is important to define what we call a 
low temperature: foUowingJ, we denote by Lc the average 
distance between kinks on a infinite facet: we define the 
low temperature regime as that in which Lc ^ L where 
L is the typical size of our island, large facets are then 
visible on the island. It was shown that L^ — §6X^( 2^^^ ) 
where a is the lattice spacing. 

The behavior of the relaxation time as a function of the 
temperature and N, the number of particles of the island, 
can be summed up with two equations corresponding to 
the high and low temperature regimes: 



the boundary of the island. Thus the germ grows, and 
eventually completes a new row on the facet. 



t: 



HT 

relaxation 



.LT 
^relaxation 



ex e^p[3E/kbT]N^ for TV > L^ (2) 

ex exp[4^/fc6r]iV for iV < L^ (3) 

Replacing the temperature dependent factors by a func- 
tion of Nc the crossover island size (where N^, = L^ <x 
exp{E / kbT)) , these two laws can be expressed as a unique 
scaling function depending on the rescaled number of par- 
ticles N/Nc-. 
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So that the relaxation timcQ is a simple monotonous func- 
tion of N/Nc, and the temperature dependence is con- 
tained in Nc- 

We will now focus on the precise microscopic description 
of the limiting step for relaxation in the low temperature 
regime. 



II. DESCRIPTION OF THE LIMITING PROCESS 
AT LOW TEMPERATURE 

A. Qualitative description 

During relaxation at low temperature, islands are 
mostly in fully faceted configurations. Let us, for in- 
stance, consider an island in a simple configuration given 
by fig. |l|. When L is larger than I, the island is not in its 
equilibrium shape (which should be more or less a reg- 
ular hexagon). To reach the equilibrium shape, matter 
has to flow from the "tips" of the island (facets of length 
I in this case) to the large facets L. In this low tempera- 
ture regime there are very few mobile atoms at any given 
time, therefore this mass transfer must be done step by 
step: the initial step being the nucleation of a "germ" of 
two bound atoms on a facet of length L and then, the 
growth of this germ up to a size L — 1 due to the ar- 
rival of particles emitted from the kinks and corners of 




FIG. 1. Configuration of the island we will consider in the 
qualitative description as well as in the quantitative one. 

This simple picture still leaves a basic question unan- 
swered: the relatively faster formation of a new row on 
a small facet would lead the island further away from its 
equilibrium shape, and yet, we observe that this never 
happens. Indeed, sometimes a germ appears on a small 
facet but it eventually disappears afterwards, whereas the 
appearance of a germ on a large facet frequently leads to 
the formation of a new row, taking the island closer to 
its equilibrium shape. 

These observations are at the root of irreversible nature 
of the relaxation, germs only grow and become stable on 
the large facet, so the island can only evolve to a shape 
closer to equilibrium. Yet, there is clearly no local drive 
for growth on large facets nor any mechanism inhibit- 
ing growth on small ones. In order to explain how this 
irreversibility comes about, we propose the following de- 
tailed description of the mechanism of nucleation and of 
growth of a germ. 

First, to create a germ, 2 atoms emitted from the cor- 
ners of the island have to encounter on a facet. The 
activation energy required for this event is obviously in- 
dependent of whether it occurs on a large facet or on a 
small facet. Once there is a germ of 2 atoms on a facet, 
the total energy of the island does not change when a par- 
ticle is transfered from a kink to the germ (3 bonds are 
broken, and 3 are created) see Fig. 0. Clearly the same 
is true if a particle from the germ is transfered to its site 
of emission or any other kink. Thus, germs can grow or 
decrease randomly without energy variations driving the 
process. An exception to this occurs if the particle that 
reaches the germ is the last one of a row on a facet; in 
that case, the energy of the system decreases by 1 bind- 
ing energy E. The island is then in a configuration from 
which it is extremely improbable to return to the previ- 
ous configuration. For this to occur, a new germ would 
have to nucleate (and grow) on the original facet. This 
event is almost impossible in the presence of the kinks 
of the first growing germ, which act as traps for mobile 
atoms. Thus, when a germ nucleates on a facet, it can 
grow or decrease without changing the energy of the is- 
land except if a complete row on a facet disappears, in 
which case it "stabilizes" . 





FIG. 2. The total energy of the island does not change 
when one particle leaves a kink on the small facet to go in the 
germ on the large facet. 

The scenario above explains why no new rows appear 
on small facets: when a germ grows on a small facet, 
since atoms come either from a small or a large facet, no 
complete row of a facet can disappear during the germ's 
growth, and thus, the island never decreases its energy. 
On the other hand when the germ grows on a large facet, 
the germ might grow or decrease, but if the size of the 
germ reaches the size of the small facet, the energy of the 
system will decrease and the system has almost no chance 
to go back its previous shape. We believe that this is the 
microscopic origin of irreversibility in the relaxation of 
this system. It should be stressed that this scenario for 
the growth and stabilization of germs is different from 
usual nucleation theory, where the germ has to overcome 
a free energy barrierQ to become stable. 

This microscopic description also shows that the lim- 
iting step for this "row by row" relaxation mechanism 
is actually the formation on a large facet, of a germ of 
the size of the small facet. This fact allows us to esti- 
mate the duration of the limiting step at each stage of 
the relaxation process. 



B. A Quantitative description 

Based on our description of the process, we propose a 
scheme for calculating the time required to form a sta- 
ble germ, i.e. a germ of size I, on a facet of size L. As 
mentioned above, the appearance of this stable germ is 
the limiting step for the formation of a new row on that 
facet. 

The idea is to describe the growth of the germ as a suc- 
cession of different island states, and calculate the prob- 
ability and the time to go from one state to another in 
terms of the actual diffusive processes occurring on the 
island surface. These states form a Markov chain, the fu- 
ture evolution of the system being essentially determined 
by the state of the system, independently of its previous 
behavior. 

As a further simplification, we consider a simple fully 
faceted island in a elongated hexagonal shape whose 
facets are of length L and I, see Fig. |^; moreover, we 
normalize every length by the lattice spacing a. 

The different states we consider are (See Fig. y also) : 
- state : there is no particle on the facets. 



- state 1 : one particle is on one end of a facet L. 

- state 2 : 2 particles are on the facet L: one of them is 
on one end of the facet, and the other one has diffused 
from an end. 

- state 3 : 2 particles are on the facet L but they are 
bonded together. 

- state 4 : 3 bonded particles are on the facet L. 

- state n : n — 1 bonded particles are on the facet L. 
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FIG. 3. Outline of the states considered in the Markov 
chain. 

The goal of this calculation is to estimate the time to 
go from state to state I + 1. Wc treat the problem as a 
discrete time Markov chain. The unit time being T2, the 
typical time for a particle with 2 neighbors to move. This 
time is in fact the smallest relevant time of the system so 
that the operation of discretization does not affect the re- 
sults. In the following, the time r^ is the average time for 
a particle with i neighbors to move: t; = i^q e'i^ . For 
clarity, we will use the term time for the discrete time of 
the Markov chain, and the term real time for the time of 
physical process. The obvious relation between the two 
time scales is : time = L^sUlini^ 

We define the parameter p — —. In the limit of small 
temperature, p is a very small quantity. Moreover, one 
can easily check that p = exp{—-^) = -p- where L^ is 
the average distance between kinks defined in the first 
section. So that the condition Lc^ L (low temperature 
regime) could be written as: ^^L ^ 1 or pi ^ 1/L < 1. 

Denote by ai the probability for the system in state i 
to stay in the state i the following step, pi the transition 
probability for system in state i to go to state i + 1, and Qi 
the transition probability for the system to go into state 
i-1. 

We have now to evaluate the different quantities ai, pi 
and qi in terms of the diffusive processes that take place 
on the island's boundary. 

We first evaluate the quantities pi,qi and ai. Let us 
assume that the states and 2 are absorbent; the aver- 



age time ni needed to leave state 1, corresponds to the 
average real time a particle stays on a facet starting on 
one of its edges, which act as traps. Since the particle 
performs a random walk, this can be readily calculated 
to be Lt2- So we should have : 



ni 



1 — ai 



= L 



(4) 



Moreover, the probability to go to state 2 is the prob- 
ability that a new particle leaves a kink and reaches the 
facet while the first particle is still on the facet, we cal- 
culate this probability in Appendix ^ where we find: 



P =1- 



2 sinh(2^(L - 1)) sinh(2^(L/2 - 1)) 



sinh(2ypL) 



sinh(2^i/2) 



(5) 



We expand expression Eq. || for small p keeping the first 
term : 



P - 2p{L - 1) + o{p) 



(6) 



We could also calculate the probability 6i, that the sys- 
tem in state 1 eventually reaches state 2 as: 



6i = aibi +pi. 



(7) 



Thus 6i 



= PI 

L 1-ai 


and using pi + a^ ^ 


P 


- 2P^^ + o(p) 


ai = 1 — 


1/L 


1 


P 1 „ L-l 



91 = 1 


we 


find: 

(8) 
(9) 


fo(p) 




(10) 



We can do the same with state 2, knowing that the 
probability for two particles to stick (state 3) is A/L, 
where A = (^^^^jf^) ti" (the calculation of this proba- 
bility is carried out in appendix ^). The time to leave 
state 2 is kL, where k is a numerical constant given by 



we obtain: 



1 l-|-4fc2 



P2 
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(see appendix |c|). Thus, 
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(11) 
(12) 

(13) 



In order to obtain a chain which can be treated analyt- 
ically, we assume that the probability to go from state 3 
to state 4, is the same as the probability to go from state 
4 to state 5 and, in general, that the probabilities to go 
from state i to state i -I- 1 are pi = Pi+i = p for i > 3. 
Similarly, we assume that ai = ai+i = a for i > 3 and 
qi = 9i+i = 9 for i > 3. 



To calculate the probabilities p, g, a, in appendix ^ 
we have calculated the average real time ip,g,a to go 
from state i to state i -f 1 assuming the average dis- 



tance between kinks and the germ is 
L/4 * T3 -I- L(L - 2)/4 * T2. Moreover 
we can calculate p, q and a : 

2(i-2 



L/2: t 



p,q,a 



Since p 



p = 2p/L - 
q = 2p/L - 



L 

4{L 



L 
2(L-2 



2)^, 



+ o{p') 
o{p') 



(14) 
(15) 
(16) 



When 2 particles are bonded on the facet (state 3), the 
probability to go to state 2 should practically be equal 
to g, we will assume this to be the case. 

So far, we have omitted the possibility that the germ 
can also nucleate on a small facet /. To take this into 
account, we have to consider new states : 

- state —1: one particle is on the facet / on one of its 
edges. 

- state —2: 2 particles are on the facet I: one of them is 
on the edge of the facet, and the other has diffused from 
an edge 

- state —3: 2 particles are on the facet / but they are 
bonded. 

- state —4: 3 bonded particles are on the facet /. 

- state —I: I — 1 bonded particles are on the facet I. 

As discussed above, if the system arrives to state —I, 
a row on a small row is completed, which is not an ab- 
sorbent state, and since this row cannot grow any further 
the system can only go back to state — Z -I- 1 or stay in 
state —I. 

The different probabilities of transition from one state 
to another in this branch of the chain are the same as the 
ones calculated before, replacing L by L Thus we have 
9-i = Pi{L => I), P-i = qi{L => I) and a_i ^ ai{L => 
I) for i > 1 except for the state —I, where we always 
have a-i — a{L ^ I), but, p_/ = 2 * q[L ^ I). In the 
following, we will use the notation : p* — Pi{L => I), the 
values of pi where we have replace L by I. 

To complete the calculation, we now have to determine 
Pq , go and ao ■ The average real time the system stays in 
state 0, assuming states 1 and —1 are absorbent, is almost 
T3/2 if we take into account that there are two kinks, one 
at each end of the facet. The probability that the germ 
nucleates on the facet L is simply j^,2r Froni this we 
deduce: 

^ "^ ;P (17) 



Pq 



ao 



90 



L + 2r 
l-2p 

L + 21' 



(18) 
(19) 



The diagram of the entire Markov chain is then given 
by Fig. ^. And we will calculate the time to go from state 



to state I + 1. 



III. ESTIMATION OF THE RELAXATION TIME 
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FIG. 4. Diagram of the entire Markov chain with the two 
branches : the upper one for the germs growing on a small 
facet and the other one for the germs growing on the large 
one, state / + 1 is the only absorbing state. 

The state / + 1 is absorbent (as discussed above, when 
the size of the germ reaches I on a large facet, the system 
cannot go back to the initial state). Let us call n, the 
average time to go from state i to state I + 1. We can 
write : 



(20) 
(21) 
(22) 

(23) 
(24) 
(25) 
(26) 
(27) 

(28) 
(29) 



The boundary condition for this process is n;+i = 0. The 
calculation of no is quite straight forward but tedious, it 
is carried out in the appendix |^. We find that in the 
limit of small temperatures, the typical real time needed 
to nucleate a germ of size I on a facet of size L is given 
by: 



with 3 < k < I 

n_2 = 1 + q2n-i + 0*211^2 +P2'^^3 
n_i = 1 + qlna + a\n^i +p*n_2 
nQ = l + qon-i + aono + poni 
ni = 1 + qino + aiUi + pin2 
"2 = 1 + 92^1 + a2"2 + P2»^3 
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= 1 + qnk-i + auk 


with 


■i<k<i^i 



r{L,l) 



rj {L + 21) 
T2 AL{L - 1) 



(__1)(/_1) + 1 



(30) 



where we have kept only the most relevant term at low 
temperatures. 



Scaling laws 

In this section we calculate the typical time required 
for an island to relax from an initial out of equilibrium 
shape. We assume that at all times, the instantaneous 
shape of the island can be characterized by the lengths 
L and I of its long and short facets respectively. Then, 
following the discussion in the previous sections, a new 
row of particles will appear on a long facet after a time 
t{L,1). Thus, calling v{L) the normal speed of the large 
facet and taking the particle size as our unit distance, we 
have: 



viL) 



1 



4t2L{L - 1) 



r(L,0 



KL + 20[(f-l)(/-l) + l] 



(31) 



The scaling properties of the relaxation time can be 
deduced by noticing that the length scales involved scale 
as N^/"^, where N is the number of atoms of the island. 
Thus we renormalize the lengths by x ^ x' = N^^''^x. 
Then, to scale out the size dependence as N grows, one 
must rescale time by i — + t' = tN^^. This is the result 
obtained ircl: at low temperature, the relaxation time is 
proportional to the number of the atoms of the island. 
But, as we expect our results for the time r equi red to 
complete a row at each stage, as given in Eq. ElO , to be 
relatively accurate, we can go beyond the scaling proper- 
ties and use it to calculate numerically the time required 
for the complete relaxation process, including the correc- 
tions arising from the lower order terms. 

In what follows, we establish the differential equations 
which permit the calculation of the full relaxation time 
of an island. As mentioned above, we still consider the 
simple island of Fig. n^, where v{L) is the normal speed 
of the facet L, and v(T) that of the facet /. We now con- 
sider L and I as continuous variables, which considerably 
simplifies the calculation. 

Conservation of the matter imposes the relation: 



Lv{L) + 2lv{l) = 



(32) 



Moreover, we can find geometric relations between 

L,l,viL),vil) : 



v{L) = V3/2 



(U 
di 



,,, /— , , Oil CLLj , 



So finally we find: 

(U 
dt 



V^t{L,1) 
dL _ 2 /L \ 1 



(33) 
(34) 

(35) 
(36) 



To integrate these equations numerically: we use 
Eq. EIQ, a nd for pl,ql and ai the exact estimation us- 
ing Eq. A9, as well as the explicit values of the quantities 
we have found: 



qi,ai and pi we have found: Eqs. |ll| - |l9[ . We start the 
integration from an island of aspect ratio of i? = 10, and 
stop it when the aspect ratio is R = 1.2. Aspect ratios 
are explicitly calculated as: 



R 



rx 




(37) 
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(38) 
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(39) 


/ / dxdy 

J J Island Surface 




(40) 



Where xq and yc give the position of the center of grav- 
ity of the island. We report in Fig. the relaxation time 
as a function of iV, the number of particles of the island 
in a log- log plot. 

We find a good quantitative agreement between the 
simulations and our predictions except at the highest 
temperatures (see below). 
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FIG. 5. Size dependence of the relaxation time for differ- 
ent temperatures in log-log bases. Filled points have been 
obtained by Monte Carlo Simulations whereas the solid lines 
have been obtained by integration of the system p3, pq. The 
agreement of the two analysis is almost perfect. 



IV. SUMMARY AND DISCUSSION 

We have considered the shape evolution of two dimen- 
sional islands as a result of the nucleation of a germ on a 
facet which then grows or decreases due to single particle 
processes. Then we recognized that the disappearance of 
a complete row is responsible for the stabilization of the 
germ, and that this is only feasible for germs growing on 



the large facets. This gives rise to an overall flux of parti- 
cles from the small facets to the large facets which leads 
the shape of the island irreversibly towards equilibrium. 
Based on this description, we have recast the formation 
of stable germs, which is the limiting step for relaxation, 
into a Markov chain in which the transition probabili- 
ties are calculated in terms of the underlying diffusive 
processes taking place on the island's boundary. Solving 
this chain yields an estimate of the time of formation of a 
new row at each stage of the relaxation. Integrating our 
results we can obtain the relaxation time for the evolution 
of an island from an aspect ratio of 10 to an aspect ratio 
of 1.2, as a function of temperature and island size. Our 
results have a rather good quantitative agreement with 
those obtained from direct simulations of the system. At 
higher temperatures, multiple nucleation processes, the 
presence of many mobile particles on the island's bound- 
ary and the failure of our hypothesis Lc^ L (Sect. II B) 
invalidate our picture and the relaxation becomes driven 
by the coarsed grained curvature of the boundary, which 
leads to MuUins classical theory. 

The description of the simple model considered in this 
work is certainly not exact and there are other effects 
that could be taken into account. Perhaps the most im- 
portant effect we have overstepped at low temperatures, 
is related to our assumption that after the germ stabi- 
lizes a single new row is formed on the large facet. The 
differential equations for the evolution of the island were 
derived from this assumption. It is clear that this is not 
correct: our estimation of the time required to stabilize a 
germ starts from a fully faceted configuration, and once 
a complete row on a small facet disappears the germ be- 
comes stable and a full row on the large facet can be 
formed. Once this row is finished, it is very unlikely that 
the island will be in a fully faceted configuration again, 
leaving at least an extra kink on the boundary. This gives 
rise to extra sources and traps for mobile particles, which 
might have affected the relaxation rate. Another issue is 
our characterization of the faceted island with only two 
facets sizes: a more detailed characterization may be rele- 
vant especially in the early stages of relaxation. It is clear 
that more accurate models for specific systems can also 
be constructed. These could take into account the depen- 
dence of the edge-diffusion coefficients on the orientation 
of the facet, as well as the dependence of emission rates on 
the local geometry. Such dependences have been studied, 
for example, bycl. In terms of the elements of description 
we use, inclusion of these effects would be achieved by 
changing the values of T2 (diffusion time) and T3 (emis- 
sion time from kinks and corners) depending on the ori- 
entation of the facets involved in each event. Thus, the 
nucleation time would depend on the facet upon which it 
happens. Such a dependence of the nucleation time may 
drive the island toward a non regular hexagon equilib- 
rium shapdl3, and would reproduce the phenomenology 
of a larger variety of materials. These changes will af- 
fect the temperature dependent prefactors in our results, 
as these depend on the temperature through T2 and T3. 



However, the size dependence of the nucleation time and 
of the relaxation time, which is where the departure from 
MuUin's theory is evidenced, would stay the same. 
From a more general point of view, only the diffusion of 
particles along the perimeter of islands has been taken 
into account in this work. In real systems, other mech- 
anisms can contribute to the transport of the matter 
which leads to relaxation: volume diffusion and transport 
through the two dimensional gas of particle surrounding 
the island. Volume diffusion is usually a much slower 
process than the other two, and can usually be neglected 
safely. On the other hand, it is well known that edge dif- 
fusion is more efficient for short trips whereas transport 
through the 2D gas is faster on long distances. Following 
Pimpinelli and Villain ina(p.l32), a characteristic length 
ri beyond which edge diffusion is less efhcient than trans- 
port through the gas can be evaluated : ri « \/IJ^t\, 
where Ds is the edge diffusion coefficient and 1/t„ the 
probability per unit time a given particle leaves the is- 
land. So that our assumptions should be valid for islands 
with a number N of particles such that : N <C DgTy. 
Since the activation energy for edge diffusion is smaller 
than the activation energy for evaporation, ri is a de- 
creasing quantity with temperature, and we expect ri 
to be very large at low temperature. Thus, this mech- 
anism is essentially irrelevant in the description of the 
evolution of nanometer structures atpLow temperatures. 
Moreover, recent experimental resultald have shown that 
supported Ag two dimensional islands relax via atomic 
diffusion on the island perimeter, without significant con- 
tribution from exchange with the two dimensional gas. 
Finally, our results are to be compared with a recent 
theoretical studyli3 concerning the relaxation of three di- 
mensional crystallites. This study also points out two 
relaxation regimes as a function of temperature. At high 
temperature the relaxation scales in accordance with the 
results derived from MuUins' theory, whereas at low tem- 
perature the relaxation time becomes an exponentional 
function of the size of the crystallites. So that the effects 
of lowering the temperature are qualitatively different for 
two dimensional and three dimensional crystallites : in 
two dimensions, lowering the temperature decreases the 
strength of the dependence of the relaxation time as a 
function of the size of the crystallites (as it crosses over 
from a iV^ dependence to a iV dependence), whereas it 
increases this the strength in three dimensions. In both 
cases, the limiting step is the nucleation of a germ on 
a facet : a unidimensional germ in two dimensions, and 
a two dimensional germ in three dimensions. The dif- 
ference stems from the fact that in the two dimensional 
case, the activation energy for the creation of the germ 
does not depend on the size of the island, it is always con- 
stant : AE, and it stabilizes when a row on a small facet 
has been removed. In the three dimensional case, this 
activation energy depend on the size of the crystallite. 
The transfer of a particle from a tip of the crystallite to 
the germ has a gain in volume energy (depending on the 
size of the islands) and a loss in edge energy of the germ 



(depending on the size of the germ) . Summing these two 
terms, an energy barrier proportional to the size of the 
crystallite appears for the creation of a stable germ. The 
exponential behavior of the relaxation time as a function 
of N is a consequence of this energy barrier dependence. 
Finally, we believe that the overall picture presented here, 
while still oversimplified, seems to be complete enough to 
provide a general picture of the processes leading to the 
shape relaxation of two dimensional islands at low tem- 
peratures. 
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APPENDIX A: CALCULATION OF THE 

PROBABILITY TO HAVE 2 PARTICLES ON THE 

FACET 

We calculate the probability P of having 2 particles on 
a facet with absorbing boundaries, knowing that at time 
t = 0, one particle is on one edge of the facet (abscissa 1), 
and that the other particle can appear on the facet with a 
probability per unit time 1/r. Relatively to our Markov 
chain, this probability is the probability the system in 
state 1 eventually reaches state 2. 

We denote by S(x, t) the probability that a particle 
on the facet at position x at time t = is still on the 
facet at time t. Then, S{x, t) satisfies the usual diffusion 
equation: 



dS{x,t) _ d^S{x,t) 



dt 



dx"^ 



(AI) 



This equation is to be solved with the conditions: 
S{x, 0) = 1 for every x s]0, L\ (i.e. we are sure to find the 
particle on the facet at time t — 0), S'(0, t) = S{L, t) ~ 
for every t (i.e. the boundaries sid es o f the facet are 
absorbent). The solution of Equation |Al| is 

S{x,t) = Y, sm^ -^e P . 

^ (2n+ l)7r L 



ri=0 



(A2) 



To take into account the appearance of particles on 
the facet, we assume the process to be Poissonian so that 
the probability to have a particle appearing at time t is 

T 

Thus the probability P that two particles are on the 
facet is: 



P 



1 



e ^ S{x,t)dt 



(A3) 



To take into account that particles can appear on the 
facet on both ends, we take r = r3/2, which holds at low 
temperatures. This leads to the expression : 



2 +°° 



P = 



TT ^ 2n + 1 

ra=0 



sin((2n+l)x) 
(52 + (2n + 1)2 



where: 



2L2 



Using formulali^ 



+ 00 



c(5,x) = E 



Dn^T3 



^ = T 



cos(fcx) 
,52 + A:2 



(A4) 

(A5) 
(A6) 



k=l 

TT cosh((5(7r — x)) 



we have 



25 sinh(<57r) 

1 sin((2n + l)x) 
'^ 2n + 1 (52 + (2n + 1)2 

n— ^ ^ 



1 

2^2 



- KH (A7) 



+ 00 

E 



C(^,x')-4CW4,2x') 



dx' 



(A8) 



So that, using Eqs. A^, A6, we finally find the following 
expression for P: 



P= 1 



2 sinh(2^(L - 1)) sinh(2^(L/2 - 1)) 



sinh(2^L) 



sinh(2^L/2) 



(A9) 



where p — —, and we have taken x = 1 (the initial 
particle on the face is at position 1 at time i = 0). We 
can expand expression Eq. A9 for small p keeping the 
first term : 



P = 2p{L - 1) + o{p) 



(AlO) 



APPENDIX B: PROBABILITY THAT TWO 
PARTICLES STICK ON THE FACE 

In this part we evaluate the probability Pa that two 
particles stick on a facet with absorbing boundaries, 
knowing that at time t = 0, one particle is on one end of 
the facet, and the other is at a position xo on the face. 
This problem can be mapped to a 2-d problem in which, 
at time t, the first particle is at position y, and the sec- 
ond particle is at position x. This virtual particle moves 
diffusively in a square of side L; starting from position 
(xo,a) (a is the lattice spacing). The quantity we are 
looking for is the probability for this virtual particle to 



reach the diagonal y = x oi the square. Thus we can 
consider the motion of the virtual particle in the triangle 
< a; < 2/ < i. We call D{x,y) the probability that a 
virtual particle starting at time t = from (x, y) leaves 
the triangle by the diagonal, V{x, y), the probability that 
this particle leaves the triangle by its vertical side, and 
H{x,y) the probability that this particle leaves the tri- 
angle by its horizontal side. We use here a continuous 
description: the discrete problem being far too difficult. 
It can be easily seen that D,V and H satisfy Laplacian 
equations: 

AD{x, y) = AV{x, y) = AHix, y) ^ (Bl) 

With the conditions : 



D{x, x) = 1 D{x, 0) = £'(0, y) = 
V{x,x)^0 F(a;,0) = V{L,y)^l 
H{x, x)^0 H{x, 0) = 1 H{0, y) = 

for Va; e [0, L] and Vy G [0, L] 

And moreover we should have : 

Dix,y) + V{x,y)+H{x,y)^l, 



(B2) 
(B3) 
(B4) 



(B5) 



which states that the particle is sure to leave the trian- 
gle since all sides are absorbent. Instead of calculating 
directly D{x,y), we will calculate V{x,y) and H{x,y). 

We first calculate Vn(x,y) and Ha{x,y) which are the 
probability that one brownian particle in a square with 
absorbing sides, respectively leaves the square by the ver- 
tical side X = L and by the horizontal side y = 0. So we 
have : 



AVnix,y)=0 AHu{x,y) = 0, 

with boundary conditions : 

Vn (x, L) = Vn (x, 0) = 

Vu{Q,y) = Vu{L,y)^l 

Ha{L.,y)^0 iJo(0,y) = 

Ha{x,Q)^l Ho{x,L) = 

for Vx (E [0, L] and Vy G [0, L] 



(B6) 



(B7) 

(B8) 

(B9) 

(BIO) 



The solution of equations Eq. B6, with these conditions 
is: 



Hn{x,y) = 



Va (x, y) 



4 Y^ sm 



{2m+l)TTX 



sinh 



(2m+l)-n-(L-y) 



TT 






m=0 

4 v^ sm 

TT 



2m +1 

{2m-\-l)TTy 



sinh(2m + l)7r 
sinhi2!2+l)££ 



m=0 



2m + 1 sinh(2m + l)n 



(Bll) 



(B12) 



One now can deduce the values of V{x, y) and i?(x, y) 
of our initial problem with a superposition of solutions 
to impose the condition : V{x,x) = and H{x,x) = : 



V{x, y) = Va (x, y) - Va (y, x) 
iJ(x, y) = Ha (x, y) - Ha (y, x) 



(B13) 
(B14) 



Using Eq. ^, ^l|, ^l|, ^Tl|, |bT|, we finally find an 
expression for D{xo, a) '■ 



D{xo,a) 



7n—0 
m—O 



(2rri+l),ra 
L 

2m+l 

(2m+l),rxo 



(2,7.+ l),rxo 



— sinh 



(2m + l),r(L-xo) 



sinh(2rM+l)7r 



2m+l 



(2Tn+l)j(L-a) 



— sinh - 



sinh(2rM+l)7r 



We now have to calculate the probability that the sec- 
ond particle is at position xq when the other one appears 
on the facet. We denote this probability P{xo)- Then 
the probability we are looking for is simply : 



A 



D{xo,a)P{xo)dxo 



(B15) 



We are actually able to calculate the exact probability 
P{xo), but then we are not able to find a simple expres- 
sion for Pa, so we prefer to make the following approx- 
imation: in the limit of small temperature, the typical 
time needed for a particle to appear on the facet is about 
T3, which is long compared to the typical time a particle 
lasts on the facet (about Lt2 for a particle that starts near 
the edge). Thus, to a good approximation, the probabil- 
ity P{xo) has reached its stationary value. Taking into 
account only the first term of the series, we have: 



Pixo) 



n , TTXo 



(B16) 



With this expression Eq. B15| becomes easy to calculate 
and we find : 



sinh ^^^^ - sinh ^ 

Pa = 1 ^ ^ 

smhTT 

/coshvr-f 1\ na n^a'^ 

^ V sinhTT / T ^ 2L2 



(B17) 
-o{l/L^) (B18) 



So that we find that the leading term of Pa is propor- 
tional to 1/L. We write this as: 



P 



A 



with A 



^( 



A-i+o(l/L') 



cosh-/ 



1 



\ sinh TT 



(B19) 
(B20) 



APPENDIX C: CALCULATION OF THE TIME 
TO LEAVE STATE 2 

We calculate in this part the average time < t > for 
two particles on the facet either to bond, or either one of 
them to reach the boundary. This time corresponds to 
the time the system stays in state 2 in the Markov chain. 
Using the description of appendix H with the virtual par- 
ticle in the triangle, we are looking for the average time 
this particle needs to leave the triangle. 



Let us call S{xo, j/o, t) the survival probability: i.e. the 
probability a particle starting at (a;o,2/o) at time t = 
is still in the triangle at time t. The average time < 
■'"(a^o, yo) > the particle stays in the triangle is given by : 



< T(a;o,yo) >= 



S{xo,yo,t)dt 



(CI) 



Moreover calling P/\{x,y,xo,yQ,t) the probability that 
one particle starting at (xo, J/o) at time i = is at position 
(a;, y) at time i, we have : 



S{xQ,yQ,t)^ 11 PA{x,y,xo,yo,t)dxdy (C2) 

Finally, we are only interested in < t{xo, 1) > since we 
know that our virtual particle starts at (xo,a), so we 
have to average the time < t{xq, 1) > over all values xq. 
To do th is we use the approximate distribution given in 
Eq. ^T^: 



<T>^ I P{xq) < t{xo, 1) > dxo 



(C3) 



So finally, using Eq. |C1|, |C2| , |C3|, we need to evaluate: 

P{xq)P/\{x, y,XQ, 1, t)dxdy dt dxQ 



<T >-- 



A 



(C4) 



We now have to calculate P/\[x, y, xq, 1, t). As before, 
we find the solution in a square, and then by superposi- 
tion, we deduce the solution in the triangle : 

PA{x,y,xo,l,t) = Pn{x,y,xo,l,t) - Pa{x,y,l,xo,t) (C5) 
Pn{x,y,XQ,l,t) = 



h E™,„=i sin ^^ sin ^ sin isp sin I^e 



(C6) 



We find < r > integrating Eq. C4 using Eq. C5, C6 over 
Xq first, and then over x and y, and finally over i, we find 



< T > = 



_ P2 ^ 1 . 2fc^ 

>- ]5^ Z^ 1 + 4fc2 ^i'^ — 

fe=i 



1 



k k{2k + l) 



(C7) 

We interested in the leading term as i ^ 00, let us define 
a function f{u) by : 



C30 _. 



k=l 



1 



k fc(2fc+l) 



f{u) is a normally convergent series, so that : 



df{u) 
du 



C30 -. 

E- -— r cos(ku 
1 + 4A;2 ^ 
fc=i 



(2A: + 1) 



fc=i 



4/c2 



1 



(2A:-f 1) 



(C8) 

(C9) 
(CIO) 



which is also convergent. The second term of the de- 
velopment is of order v^ or sma ller, so that integrating 
Eq. |C1C, and using it in Eq. C7, we find: 



< T >= 



2i 






4fc2 



2- 



1 



(2fc+l) 



0(1/L) 
(Cll) 



Thus < r > is proportional to L, and we write: 

<t>=hLt2 + 0{1/L) (C12) 



with K, 



7^2 Z^ 1 + 4fc2 



1 



(2fc + l) 



(C13) 



APPENDIX D: CALCULATION OF THE TIME A 

PARTICLE NEEDS TO GO TO CENTER OF THE 

FACET 

We now calculate the average real time tp,g,a to go 
from state i to state i + 1 assuming the average distance 
between the kinks or corners from which particles are 
emitted and the germ is L/2. This again can be posed 
as a Markov chain where : 

- state : no particle is on the facet; 

- state 1 : one particle coming from a kink is in position 
1 on the facet; 

- state 2 : the particle is at position 2; 

- state k : the particle is at position k; and state L/2 is 
absorbent. 

Again we use T2 as the unit time. Using the same 
definition for the quantities pi, qi and ai as in the main 
text, we have: 



P^ = 1/2 

q^ - 1/2 
with i > 1 



(Dl) 
(D2) 
(D3) 



To calculate po and ag we know that the real time an 
atom needs to leave a kink is ra, and as there are two 
kinks at the ends of the facet, 7^/2 is, to a good approx- 
imation, the average time to leave state 0. Thus that we 
find: 



Po = 2/3 
ao = 1 - 2p 



(D4) 

(D5) 



Calling rii the average time to go from state i to the 
absorbent state L/2, we have the equations : 



no = 1 + ao"-o +^0?^! 

rih = i -\ 

2 2 

with k > 1 
nL/2 = 



(D6) 

(D7) 

(D8) 



and ni^/2 — 0. Summing equations D7 from k ==1 to j. 



and then from j = 1 to L/2 — 1, and using Eq. D8 
finds 



L-2 



L-2 



ni 



-no 



And Eq. D6 permits to obtain ?io : 
L L{L-2) 



no 



2po 



(D9) 



(DIO) 



Going back to real time, one finds the average time r a 
particle needs to leave a kink and reach the center of the 
facet is: 



L L{L-2) 



'^p.q.a 



4 ^ 



-T2 



(Dll) 



APPENDIX E: CALCULATION OF iVo 

To carry out the calculation, we can first calculate the 



Uk for fc > 1. Noting that p = q, we have from Eq. 28 

{uk+i - rik) = {rik - rik-i) - 1/p, 



(El) 



so that summing from k=4 to j, and then from j=4 to I, 
and using nj+i = 0, we find: 

{I - 1)^3 ^{l~ 2)n2 -f 1/P(^^— ^^^— ^) (E2) 

Inserting this result in Eq. E^ we find n2 as a function of 
rii, in the same way, using Eq. Eq, one can deduce the 
equation for rii : 



ni (1 — ai — PiB) ~ 1+ piA + qiiio 



where : 



1 -U P2 1-2 

A^ ^ ' 



I- 0.2- P2JZ^ 



B 



92 



I- a2~ P2 



1-2 
l-l 



(E3) 

(E4) 

(E5) 
(E6) 



We now calculate the case fc < 0. 
Calling rrij — n^j for every j, and summing Eq. UV from 
j to z — 1 , one finds : 



rriz — ruz-i — rrij — ruj-i 



(E7) 



and Eq. |2C| yields: nii — mi-i — ^-t. Using this in Eq. E7 
and taking j = 3, one finds : 

l~3 1 
p* 2p* 



m3 ~ 1712 = 



(E8) 



10 



Using Eqs. E3, E4, one finds 



1 
n-i = — 
91 



Pi 

9192 



-4^(1 -5/2)+ no. (E9) 

9i92P 



One can know obtain the value of uq from Eq. E2 , E9 
andpS : 



Po(l- 
1+Po 



9i 



1 ~ «! — Pli? 



)"-0 



1 — Qfl 
1 



PiB 



"90 — 
.9i 



PI 



P1P2 



(l - 5/2) 



9i92 9i92P 
Then, using Eq. P-|l9|] we can calculate A and B 
L (/-!)(/ -2) 



(ElO) 



A = 



Ap[k 



m-i) + i 



jLHl-l) 



L(L-2)(i-2)(i-l) 
4 



ij:-m-i) + i 



o(l) 



B = 1 - 



(^-1)G-1) + 1 



And finally, the expression of no is: 



no = 



1 {L + 21) 
'p^ 4:L{L- 1) 

+0(l/p) 



(^-1)(^-1) 



(Ell) 



We report here only the leading term: the expressions 
of the quantities Pi,qi and at permit us to calculate no 
up to order 1, but these terms are truly ugly and it does 
not seem relevant to give them here. 

Going back to real time, we find that the time t{L, I) 
to nucleate a new row on a facet is given by : 
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